Disorder effects at a nematic quantum critical point in d-wave cuprate superconductor 
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A d-wave high temperature cuprate superconductor exhibits a nematic ordering transition at zero 
temperature. Near the quantum critical point, the coupling between gapless nodal quasiparticles 
and nematic order parameter fluctuation can result in unusual behaviors, such as extreme anisotropy 
of fermion velocities. We study the disorder effects on the nematic quantum critical behavior and 
especially on the flow of fermion velocities. The disorders that couple to nodal quasiparticles are 
divided into three types: random mass, random gauge field, and random chemical potential. A 
renormalization group analysis shows that random mass and random gauge field are both irrelevant 
and thus do not change the fixed point of extreme velocity anisotropy. However, the marginal 
interaction due to random chemical potential destroys this fixed point and makes the nematic phase 
transition unstable. 
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I. INTRODUCTION 

One important reason that high-temperature cuprate 
superconductors are hard to understand is that they have 
very complicated phase diagram. The competitions and 
transitions between different phases give rise to many 
unusual properties, and hence have attracted consider- 
able theoretical and experimental efforts in the past two 
decades. Among the widely studied competing orders, 
the various phase of the anisotropic electronic liquid are 
of particular interests. Kivelson, Fradkin, and Emery 
proposed that due to the local electronic phase separa- 
tion, a number of novel electronic liquid crystal phases 
can exist in a doped Mott insulator [l|. The simplest 
of such phases is the electronic nematic phase, in which 
the rotational symmetry is broken but the translational 
symmetry is preserved. In recent years, the nematic or- 
dering phase transition has been investigated extensively 
0, HI . The resistivity anisotropy observed by Ando et 
al. in two types of cuprate superconductors provided the 
early evidence for the predicted nematic phase [4(. More 
recently, the neutron-scattering experiments performed 
in YBa2Cu3 06.45 also pointed to the existence of nematic 
phase [5j . Further evidences came from the observed in- 
plane anisotropy of the Nernst effect in the pseudogap 
region of YBa2Cu30j, @ and from the scanning tunnel- 
ing microscopy experiments performed in the pseudogap 
region of Bi2Sr2CaCu20s+,5 0- Interestingly, there are 
also compelling experimental indications for the existence 
of nematic phase in Sr3Ru20_7 [|| and in newly discovered 
iron-based superconductor Q. 

In the language of field theory, the nematic phase has 
an Ising-type order parameter that can be represented by 
a real scalar field <f>. However, the dynamics of the system 
close to the critical point can not be fully described by 
an effective 4> 4 theory [l(| when there are itinerant elec- 
trons. The interaction between quantum fluctuation of 
the nematic order parameter and the itinerant electrons 
has to be treated carefully. Close to the quantum crit- 
ical point, this interaction becomes singular, which was 



shown to be able to produce highly unusual, non-Fermi 
liquid like, behaviors [TTl - fl4j |. The nematic physics is in- 
timately related to Pomeranchuk instability and has also 
been investigated from this point of view |15l [l6j . 

Besides the pseudogap phase of underdoped cuprates, 
it is also interesting to study the nematic transition that 
occurs in the d-wave superconducting phase [I7I - I21I ] . This 
is a new example of quantum phase transitions happening 
in the superconducting dome [22f-[24|, which is a widely 
studied topic. In the superconducting phase, the nematic 
order parameter interacts strongly with the gapless nodal 
quasiparticles, which are the low-energy excitations of a 
e?-wave superconductor. This interaction remarkably af- 
fects the dynamics of both nodal quasiparticles and ne- 
matic order parameter. An early work of Vojta et al. (l7j 
presented a detailed renormalization group (RG) analy- 
sis of various types of Yukawa couplings in the d-wave 
superconducting phase, including nematic type coupling. 
More recently, Kim et al. studied the effect of quantum 
fluctuations of nematic order parameter on the spectral 
properties of nodal quasiparticles (18(. 

In actual <i-wave cuprate superconductor, the gapless 
nodal quasiparticles have a Fermi velocity vp and a gap 
velocity va, which are not equal. Indeed, the ratio va/vf 
may be as small as 1/20 [251 . This small ratio plays an 
important role because it appears in a number of observ- 
able quantities [25}. For instance, it was found (2(| that 
the dc thermal conductivity contains the large inverse of 
this ratio as 

« k A(^ + v JL) (i) 

at nearly zero temperature. This easily accessible ma- 
terial property is universal — it is independent of the 
amount of disorder [261 ] . This universality was confirmed 
by transport measurements [27]. Since the inverse of 
va/vf is so large it completely dominates k/T. 

An interesting problem is how the velocity ratio is in- 
fluenced by the nematic phase transition. Recently, Huh 
and Sachdev [l9| studied this problem by making a care- 
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fill RG analysis within an effective field theory of nematic 
ordering transition. They found that va/vf flows to a 
fixed point with va/vf 0, i.e., the inverse velocity 
ratio vf/va diverges. Therefore, the nematic ordering 
transition in ei-wave superconductor is accompanied by 
the appearance of an extreme velocity anisotropy. Since 
the diverging velocity ratio vf /va enters various physical 
properties, the predicted extreme anisotropy should have 
observable effects. In particular, the low-temperature dc 
thermal conductivity is expected to be significantly en- 
hanced near the critical point. By using a Boltzmann 
equation approach, Fritz and Sachdev [21] calculated 
the thermal conductivity enhancement near the nematic 
quantum critical point due to the divergence of vf/va- 
If this enhancement were observed in transport experi- 
ments at certain doping concentration, this would serve 
as an important evidence for the existence of nematic 
transition in d-wave cuprate superconductor. 

When studying the low-temperature transport proper- 
ties of an interacting electron system, it is hardly possible 
to ignore the disorder effects. First of all, the fermions 
are always scattered by certain amount of disorder in any 
realistic physical system. Moreover, although the elastic 
scattering due to quenched disorder is less important at 
high temperature, it dominates over the inelastic scatter- 
ing due to inter-particle interactions at very low temper- 
ature. The disorder effects should be taken into account 
when calculating the low-temperature thermal conduc- 
tivity. The nematic order parameter fluctuation can lead 
to significant enhancement of thermal conductivity only 
when the fixed point of extreme velocity anisotropy is 
stable against disorder scattering. If the fixed point is 
changed or even destroyed by disorder, the thermal con- 
ductivity enhancement will not occur in practice. It is 
therefore crucial to examine the disorder effects on the 
RG flow of fermion velocities, especially on the stability 
of extreme anisotropy of velocities. 

In general, the disorders coupled to gapless nodal 
quasiparticles in d-wave superconductor can be divided 
into three types: random chemical potential, random 
gauge field, and random mass. The difference comes from 
the different Pauli matrices used to define the fermion- 
disorder interacting terms. The effects of these disorders 
on the low-temperature transport properties of nodal 
quasiparticles have been discussed extensively [H, [29| . 
These disorders will alter the RG flows of fermion ve- 
locities. On the other hand, the RG flow of strength 
parameters of fermion-disorder couplings are determined 
by fermion velocities, and thus should be calculated self- 
consistently with the flow of fermion velocities. 

In this paper, we present a RG analysis of the inter- 
play between nematic order parameter fluctuation and 
disorder scattering. We derive a series of coupled RG 
equations of fermion Fermi velocity Dp, gap velocity va, 
and disorder strength parameter g. In the cases of ran- 
dom mass and random gauge field, the corresponding 
disorder strength parameters both flow to zero at low 
energy. Therefore, these two kinds of disorders do not 



change the flow of fermion velocities and hence the fixed 
point of extreme velocity anisotropy is stable. However, 
the strength parameter of random chemical potential re- 
mains a constant even down to the lowest energy, and 
thus is able to modify the fermion velocities significantly. 
We found that the fermion velocities do not flow to any 
fixed points, but indeed oscillate rapidly between posi- 
tive and unphysical negative values. This implies that 
the extreme anisotropy fixed point is destroyed and the 
nematic phase transition may become unstable due to 
random chemical potential. 

In Sec. UH we define the model action in the presence 
of both nematic order parameter fluctuation and disor- 
der. In Sec. IIIH we calculate the fermion self-energy 
corrections due to nematic order parameter and disorder 
scattering. The fermion-disorder vertex corrections due 
to nematic and disorder interactions are also computed 
in this section. In Sec. IIV1 we make the RG analysis and 
obtain the self-consistent RG equations for fermion veloc- 
ities and disorder strength parameter. These equations 
are solved both analytically and numerically in Sec. [V] 
From the solutions, we found that the fixed point of ex- 
treme velocity anisotropy is not changed by random mass 
and random gauge field. However, the random chemical 
potential destroys this fixed point and indeed makes the 
nematic phase transition unstable. In Sec. I VII we briefly 
summarize the results obtained in this paper and discuss 
the possible experimental detection of the predicted ex- 
treme velocity anisotropy. 



II. MODEL 



We start from the following action 



S = S, 



Srt 



(2) 



where the free action for nodal quasiparticles is 
G? 2 k duj 



S 



(2tt) 2 2tt 
d 2 k duj 
(2tt) 2 2tt 



ip\ a (-iu + v F k x r z + v A k y T x )^ la 



V2a(- iuj + V F k y T Z + VAk x T X )lp 2 a, (3) 



where t^ x ' v ' z ^ denote Pauli matrices. The linear disper- 
sion of Dirac fermions originates from the d^2_ y 2-wave 
symmetry of the energy gap of cuprate superconductor. 
Here, the spinor ip\ represents nodal quasiparticles ex- 
cited from the (§, ■§) and ( — S, — |) nodal points, and i/4 
the other two nodal points [Tg . The repeated spin index 
a is summed from 1 to Nt, the number of fermion spin 
components. The ratio va/vf ~ 1/20 between Fermi 
velocity and gap velocity is determined by experiments 
[25} . The effective action S<p describes the Ising type ne- 
matic order parameter, which is expanded (for notational 
simplicity) in real space as 

* = / d^dr{\(dr^f + ^(V0) 2 + ^ 2 + g0 4 },(4) 
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where t is imaginary time and c is velocity. The mass 
parameter r tunes the nematic phase transition with r — 
defining the quantum critical point. The parameter 
uq is the quartic self-interaction strength. The nematic 
order parameter couples to nodal quasiparticles via the 
Yukawa term 

V = / d 2 xdT{\ Q ^\ a T^ la + 4 a T X ^ 2a )}. (5) 



Following Huh and Sachdev fl9| . we now perform the 
RG analysis in the framework of a 1/Nf expansion. The 
inverse of the free propagator of the nematic order pa- 



rameter field behaves as q 



After taking into account 



the polarization effects, there will be an additional lin- 
ear g-term. At low energy regime, the g-term dominates 
over the g 2 -term, which then can be neglected. Near the 
quantum critical point, we keep only the mass term and 
assume that <f> — > </>/Ao and r — > NfrX^, leading to 



s = s. 



d xdr- 



After integrating out fermion degrees of freedom, the ef- 
fective action for the scalar field becomes 



S,p 



d 3 q 
(2tt) 3 



[r + n(g)Mg)| 2 +0(0 4 ). 



(7) 



The lowest-order Feynman diagram for the polarization 
function is shown in Fig. [T] and symbolizes the integral 



n(q,e) 



d 2 k duj 



(2tt) 2 2tt" l ' ^ 
where the free fermion propagator is 

1 



Tr[T*G°(k, W )r*G°(k + q,u; + e)], 



-iuj + VFk x T z + v\k v T x 



(8) 



As shown previously [l9j , the propagator for the nematic 
order parameter is given by 



n(q,e) 
1 



(e 2 + v 2 F q 2 x ) 



Wvfva (e 2 
1 



(c 2 



16vfva (e 2 



(9) 



in the vicinity of nematic quantum critical point r = 0. 

Disorders are present in almost all realistic condensed 
matter systems and play important roles in determining 
the low-temperature behaviors. In the present problem, 
the nodal quasiparticles can interact with three types 
of random potentials, which represent different disorder 
scattering processes. According to the coupling between 
nodal quasiparticle and disorders, there are three types of 
random fields in <i-wave superconductors: random mass, 
random chemical potential, and random gauge potential. 
All these types of disorders have been investigated in the 




FIG. 1: The polarization function for nematic order param- 
eter. The solid line represents the fermion propagator and 
wavy line represents the boson propagator. 




(a) 



(b) 



FIG. 2: One loop fermion self-energy correction due to (a) 
nematic order parameter fluctuation and (b) disorder. The 
dashed line represents disorder scattering. 



contexts of <i-wave cuprate superconductor [28|,|29(, quan- 
tum Hall effect '■'>{) , and graphene [HI, HH . In the gen- 
eral analysis to follow, we shall consider the three types 
of disorders. 

The fermion field couples to a random field A(x) as 



d 2 xV' t (x)r^(x)A(x) 



(10) 



For random chemical potential, the matrix T is T = I. 
For a random mass it is, T — t v , and for a random gauge 
field r = (t x , r 2 ). The random potential ^4(x) is assumed 
to be a quenched, Gaussian white noise field with the 
correlation functions 

(A(x)}=0; (A(x 1 )A(x 2 ))= 3 « 2 5 2 (x 1 -x 2 ). (11) 

The dimensionless parameter g represents the concen- 
tration of impurity, and the parameter i>r measures the 
strength of a single impurity. It will be convenient to 
redefine the random potential as ^4(x) — > v-pA(x), and 
then write the fermion-disorder interaction term as [3l| 



S dis = vr J d 2 xV/(x)IV(x)A(x), 

with the random potential distribution 
(A(x))=0; (A(x 1 )A(x 2 ))-^ 2 (x 1 



x 2 



(12) 



(13) 



Now the RG flow of disorder strength can be calculated 
by studying the vertex correction to the fermion-disorder 
interaction term. After a Fourier transformation, the cor- 
responding action has the form 



J d 2 kd 2 kidwV t (k,o;)r^(ki, u)A(k - ki)(14) 



This action will be analyzed together with the actions (3), 
(6), and (7). In order to perform perturbative expansion, 
we assume that g and i>r are both small in magnitude, 
corresponding to the weak disorder case. 



III. FERMION SELF-ENERGY AND 
FERMION-DISORDER VERTEX CORRECTIONS 

According to the Dyson equation, the interactions in- 
duce a self-energy correction to the free propagator of 
Dirac fermion, yielding 

G^ 1 (k,w) = —iu + v F k x r z + v/\k v T x 

-E nm (k, lo) - E dis (k, w), (15) 

where self-energy functions E nm and Edi s come from ne- 
matic ordering and disorder scattering, respectively. To 
the leading order, the corresponding Feynman diagrams 
are presented in Fig. [2j 

The nematic self-energy E nm has already been ob- 
tained by Huh and Sachdev [lfj, who found that 



dE nm (k,w) 



din A 
where 

Ci = 



= Ci(-iu) + C' 2 v F k x T z + C 3 VAk v T x , (16) 



r-2-K 

dx / dO 
Jo 



Co = 



c. 



/ 

x 2 - cos 2 9 - (va/v f )' 2 sin 2 9 
X {x 2 + cos 2 9 + {va/vf) 2 sin 2 (9) 2 

cos 2 ( 



(x 2 + cos 2 9 + {va/vf) 2 sin 

Nfir 3 J 

(va/vf) 2 sin 2 



G(x,9)(l7) 



x 2 + cos 2 i 



(x 2 + cos 2 9 + {va/vf) 2 sin 2 9) 2 



g- 1 = 



x 2 + cos 2 I 



\J x 2 + cos 2 9 + {va/vf) 2 sin 2 



+ 



a; 2 + sin 2 ( 



x 2 + sin 2 + (va/vf) 2 cos 2 # 



(20) 



The computational details of C± t 2,3 are presented in the 
appendix. 

The fermion self-energy due to disorder Edi s (iw) can 
be computed as 

•/2i 



f d 2 k 

s dis (iw) = gw r y ^-^rG^(k,a;)r 



5 w 2 



-iuj In A. 



(21) 



From this expression, we know that Edis(i^) has the same 
result for all possible expressions of T. Another impor- 
tant feature is that Edi s («w) is independent of momen- 
tum, which reflects the fact that the quenched disorder 
is static. It is easy to have 



riE dis (^j) 
dlnA 



(22) 




(a) 



(b) 



FIG. 3: Fermion-disorder vertex correction due to (a) nematic 
order parameter and (b) disorder parameter. 



where 



C„ 



2ttvfva 



(23) 



The fermion-disorder interaction parameter vr is also 
subjected to RG flow. To get its flow equation, we need 
to calculate the fermion-disorder vertex corrections. For- 
mally, the vertex correction has the form 



v r T' = v r T + V nm + Vdi s , 



(24) 



where V nm represents the vertex correction due to ne- 
matic order parameter fluctuation and Vdis represents the 
vertex correction due to disorder interaction. The corre- 
sponding diagrams are shown in Fig. [3J They will be 
calculated explicitly in the following for all three kinds 
of disorders. 



1. Random chemical potential 

We first calculate the vertex correction due to nematic 
ordering. To this end, we employ the method proposed 
by Huh and Sachdev [l9j. At zero external momenta and 
frequencies, the vertex correction is expressed as 



Km = vr 



(2tt) 3 

There is an useful formula [1 

dVam 



where 



dlnA 



1 



2tt 



H(Q) = —t x 



*T^3 / dx I d9H(Q), 



1 



(25) 



(26) 



I 



Nf ( — ivFX -\- vf cos 9t z -\- va sin 9t x ^j 
1 . 1 



(—ivpx + vp cos8t z + VAsin9r x ) H(Q) 

(27) 

Here, the matrix I corresponds to the coupling between 
Dirac fermion and random chemical potential. It will be 
replaced by t v in the case of random mass and r I,z in 
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the case of random gauge field. After straightforward 
computation, we have 



rflnA 



= C 5 v r I, 



(28) 



where 



2ts 



N fK 3 J-oo Jo 

(x 2 - cos 2 9 - (va/vf) 2 sin 2 9) 
X (x 2 +cos 2 9 + (v A /v F ) 2 sin 2 9) 2 [X ' ' 
= -Ci. (29) 

The vertex correction due to averaging over disorder is 

—^IG°(o;,p)«rIGO(a;,p + k)L (30) 

Again, the matrix I should be replaced by certain Pauli 
matrix in the case of random mass or random gauge field. 
Taking the external momentum k = and keeping only 
the leading divergent term, we have 



rflnA 



= C r v r l, 



where 



v 2 g 

°r — -z — <~'g 

Zttvfva 



2. Random mass 



(31) 



(32) 



The calculation of vertex correction in the case of ran- 
dom mass parallels the process presented above, so we 
just state the final result. The nematic ordering induced 
vertex correction is 



dV, r 



d\nA 



C 6 v r T v , 



(33) 



where 



OO /»2tT 



C 6 = 2 ±^fl I dx I d9 



NfK 3 J-oo JO 

(x 2 + cos 2 9 + (v A /v F ) 2 sin 2 9) 
(x 2 + cos 2 9 + (va/vf) 2 sin 2 9) 2 
= C3 — C\ — C-x. 

The disorder induced vertex correction is 
dV dh 



G(x,9) 



dlnA 



-C r (v r T y ), 



where 



°r — -z — o g 

Zttvfva 



(34) 



(35) 



(36) 



3. Random gauge potential 

The random gauge potential has two components, 
characterized by t x and r z respectively. For the t x com- 
ponent, the nematic ordering contribution to vertex cor- 
rection is 



dV n 



dlnA 



= C iA vrT x , 



(37) 



where 



_ _2(^m rdx[ " dl) 



oo JO 



NfTT 3 

(x 2 + cos 2 9 - (va/vf) 



2 sin 2 



(x 2 + cos 2 9 + (va/vf) 2 sin 2 9) 2 



G(x,9) 



= -Ci. 
For the t z component, we have 

dVn-m 



dlnA 



C 4B v r T z , 



(38) 



(39) 



where 



C, 



2(v a /vf) 

NfTT 3 



dx / d9 

-oo JO 
2 „„„2 a i l~. /.. \2„:„2, 



(x 2 — cos 2 9 + (va/vf) 2 sin 



,G(x,9) 



(x 2 + cos 2 9 + (va/v f ) 2 sin 2 9) 2 ' 
= -C 2 . (40) 

The disorder contribution can be calculated similarly. 
To both the t x and t z component, we have 



so that 



VdisM = finite, 



dVdxs 



dlnA 

for both two components. 



(41) 
(42) 



IV. RG EQUATIONS FOR FERMION 
VELOCITIES AND DISORDER STRENGTH 

In order to perform the RG analysis of the fermion 
velocity and disorder strength, it is convenient to make 
the following scaling transformations [l9|, [H| 



(43) 
(44) 

lMk,w) = Vi, 2 (k'^')e^« (4 -^ )dZ , (45) 
0(q,e) = <j ) '(c{,e l )e^>-^ dl 1 (46) 

where i = x,y and b = e~ l with I > 0. The parameters rjf 
and rib will be determined by the self-energy and nematic- 
fermion vertex corrections. Note the energy is required to 
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scale in the same way as the momentum, so the fermion 
velocities are forced to flow under RG transformations. 

In the spirit of RG theory 33] , to specify how a field 
operator transforms when the energy and momenta are 
re-scaled, the standard method is to require that its ki- 
netic term remains invariant. In the present problem, 
however, the random potential A(x) does not have an 
own kinetic term. In order to find out its scaling behav- 
ior, we write the Gaussian white noise distribution in the 
momentum space as 



(A(k 1 )A(k 2 ))= 9 5 2 (k 1 



(47) 



When the momentum k becomes 6k, the delta function 
is rescaled to 

5 2 (ki + k 2 ) -> <5 2 (6k! + bk 2 ) = b~ 2 8 2 {k 1 + k 2 ). (48) 

If we require that the disorder distribution Eq. (|4"7| is 
invariant under scaling transformations, then the random 
potential should transform as 



A(k) -> b _1 A(k). 
Now we have to assume that 

A(k) - A'(k')e l . 



(49) 



(50) 



According to the RG technique presented in [33j |. the 
momentum shell between 6A and A will be integrated 
out, while keeping the —iu> term invariant. From the 
nematic ordering and disorder contributions to fermion 
self-energy function, we have 

f bA A A 

/ d 2 kdwifr {—ioj — C\ (—ioS) In — + C g (— iu>) In tj^P 



d 2 kduip f (-iu;)[l + (C g - d)^ 
d 2 kduj^ (-iw)e (c ' 9 - c ' l) V- 



(51) 



After the scaling transformation, this term should go 
back to the free form, so that 



Vf=C g -C 1 . 



(52) 



The kinetic terms should also be kept invariant under 
scaling transformation, which leads to 



= (Ci-C a -C>F, 
^ = (d-Ca-C^A. 



(53) 
(54) 



Based on these expressions, the ratio between gap veloc- 
ity and Fermi velocity is 



d(v A /v F ) 
dl 



(C 2 -C 3 )(v A /v F ). 



(55) 



The disorder strength g appears in the above expres- 
sions. Due to the interplay of nematic ordering and disor- 
der, this parameter also flows under RG transformation. 
The flow equation depends on the type of disorder, which 
will be studied in the following. 

We first consider the case of random chemical poten- 
tial. The bare fcrmion-disorder action is 



J d 2 kd 2 k 1 du)^(k,u})Tilj(k 1 ,uj)A(k-k 1 ). (56) 

Including corrections due to nematic and disorder inter- 
actions yields 

r&A 



A 
bA 



/OA 
d 2 kd 2 kidu^(\s.,u)[vrl - C\v r I\n 

+C g v r Un^(k 1 ,iu)A(k-k 1 ) 

/ d 2 kd 2 k 1 dw^ t (k,a;)wrl[l + (G g - Ci)J] 
x^(ki,«)A(k-ki) 
/ d 2 kd 2 kidwV t (k,w)w r Ie (c ' 9_c ' l) ' 



xV(ki,w)A(k-ki). 



(57) 



After redefining energy, momentum, and field operators, 
we have 



d 2 k'd 2 k'idu>'i) n (k',L)')vrl 

xe( c °- ClS > l i)'{k! u w')e- r >t l A'{k! - k[). (58) 

Since rjf = C g — Ci, it is easy to obtain the following RG 
flow equation for vr, 



dvr 

~dT 



0. 



(59) 



Apparently, the parameter vr does not flow and thus can 
be simply taken to be a constant. 

In the case of random mass, the flow equations for 
fermion velocities have the same expressions as Eq.(53) 
and Eq.(54). However, the flow equation for disorder 
strength is different from Eq.(58), and has the form 

^ = (C 3 - C 2 - 2C g )v r , (60) 

which couples self-consistently to flow equations of 
fermion velocities. 

Following the steps presented above, we find the follow- 
ing RG equations in the case of random gauge potential 



= (Ci-C 2 -C gi )v F , 
^ = (d - C 3 - C gi )v A , 



(61) 
(62) 



which couple to the flow equations of disorder strength 
dvn 



dl 

dVT2 

~dT 



[(d - C gl ) - C 8 ]vn, (63) 
[(Ci - C g2 ) - C 2 }vr2, (64) 



7 




0.72 



ID"' 10"™ 10" 50 10" 30 4 10-'° 



1.0 

, 0.8 
V IV 

r ro 

0.6 
(b) 0.4 
0.2 




v Iv „ 
r ro 



1.0 

0.8 

3 

0.6 
(C) 0.4 
0.2 



io" 90 10' 7 



-g=iO" 7 
-g=i0' 5 

g=i0' 3 

-g=10" 2 



10' 



10' 



10 



10 



10 



10" 



10 4 10 

e 



FIG. 4: (a) vr for random mass; (b) vr for random gauge 
potential t x component; (c) vr for random gauge potential 
r z component. 



where 



c <9 { = 12 
9z 2wvfVa' 



(65) 



Here, the equations denoted by i = 1,2 correspond to 
the t x and r 2 components of random gauge potential, 
respectively. 



V. STABILITY OF EXTREME ANISOTROPY 
AGAINST DISORDERS 

The RG flows of fermion velocities Dp and va with 
growing scale / can be obtained by numerically solving 
the corresponding coupled equations with the initial val- 
ues Vfo, vao, an d vro- First of all, in the clean limit 
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FIG. 5: The flows of vf (upper one) and va (lower one) in 
the case of random mass. 



g = 0, the equations reduce to that obtained by Huh and 
Sachdev. In this case, it was already known that an ex- 
treme anisotropy of fermion velocities, i.e., va/vf 0, is 
caused by nematic order parameter fluctuation. The ef- 
fects of various disorders on this fixed point will be trans- 
parent when the dimensionless parameter g is increased 
smoothly. 

We first consider the case of random mass. As I 
grows, vr first increases and then decreases, eventually 
approaching zero, as shown in Fig. @|a). Although the 
RG equations for fermion velocities are modified by scat- 
tering due to random mass, the disorder parameter vr 
flows to zero as I — > oo. Apparently, the random mass is 
irrelevant in the present problem. As can be easily seen 
from see Fig. the fermion velocities vf and va both 
decrease as I grows. More concretely, va goes down to 
zero rapidly, but vf decreases much more slowly and ac- 
tually approaches a finite value. These results imply the 
existence of extreme velocity anisotropy with va/vf 
in the presence of random mass. 

We next discuss the case of random gauge potential. 
The flows of disorder strength vr with growing / are 
shown in Fig. EJb) for component t x and in Fig. Htc) for 
component t z . For both components, the correspond- 
ing vr decrease as I grows and finally approaches zero as 
/ — y oo. Similar to the case of random mass, the random 
gauge potential makes no important contributions to the 
flow of fermion velocities. Therefore, as in the clean limit, 
both vf and i>a decreases with I until approaching zero 
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FIG. 6: The running vf and va for the t x component of 
random gauge potential. The running of fermion velocities 
for component t z are very similar to this case, and thus are 
not shown. 



for component t x depicted in Fig. [6] the flows of vf and 
va for component r z with growing I will not be shown 
since they are similar to those in the case of component 
t x . It is obvious that random gauge potential does not 
change the extreme velocity anisotropy. 

Unlike random mass and random gauge potential, the 
disorder strength parameter vy in the case of random 
chemical potential does not flow with / and thus should 
be kept as a constant. As such, the influence of scat- 
tering due to random chemical potential can not be ne- 
glected and indeed the flows of velocities vf and va de- 
pend heavily on the magnitude of vr- At first glance, 
the flow equation of velocity ratio va/vf is independent 
of disorder strength vr, as shown in Eq. (55), and thus 
appears to have fixed point at va/vf — as in the clean 
limit. However, this solution is artificial. In the present 
problem, the flow equation of va/vf is derived from the 
more fundamental equations of va and vf, and therefore 
is reliable only when va and vf have well-defined fixed 
points. If the equations for va and vf have no fixed 
points, the equation of va/vf becomes meaningless. 

To see the effect of random chemical potential, we 
make a qualitative analysis based on the RG equations 
of fermion velocities va and vf ■ The fixed points can be 
obtained by requiring that 



dvp 
If 



(Ci - C 2 )v F 



2ttva 



0, 



(66) 



(Ci - C 3 )v A 



Vr9 
2ttvf 



= 0. 



(67) 



We assume that v* F and v* A correspond to the fixed points. 
If both v* F and v* A are finite, then the above equations 
imply that (C\ — C2)v* A v* F = (C\ — C 3 )v* A v* F , which can 
not be satisfied since v* A ^ 0. If v A = 0, then 



27rtft(Ci-<7 2 ) - 



(68) 



From the expressions for C\ and C2, this implies that 
I cx l/(v A ) 2 , which is clearly inconsistent with the as- 
sumption of v A — 0. Before going to v* F = case, we 
define Di = (vf / fA)Ci, % = (1,2,3), then the new forms 
of Eq. (IBBl and Eq. (|S71) become 



{Dx - D 2 )va 



v F 



Vr9 

2lTVA 

via 

2-KVp 



(69) 
(70) 



If v* F = 0, by both analytical and numerical analysis, we 
found that these equations have no solution. 

In summary, as shown by the above analysis, there is 
no fixed point of the fermion velocities vf and va when 
the fermions interact with random chemical potential. 
Therefore, vp and va do not approach any stable val- 
ues at the low energy regime. Alternatively, straight- 
forward numerical calculations show that they oscillate 
rapidly between positive and unphysical negative values 
as I grows. In this case, the extreme velocity anisotropy 
fixed point is destroyed. We interpret the occurrence of 
unphysical negative velocities as a signature of the in- 
stability of nematic phase transition in the presence of 
random chemical potential. 



VI. SUMMARY AND DISCUSSION 

In summary, we have examined the disorder effect near 
the critical point of nematic phase transition in d-wave 
cuprate superconductor. We considered three types of 
quenched disorders that couple directly to the gapless 
nodal quasiparticles: random mass, random gauge field, 
and random chemical potential. By means of a RG anal- 
ysis, we have derived a series of self-consistent flow equa- 
tions for fermion velocities and disorder strength. It was 
found that the fixed point of extreme velocity anisotropy 
due to critical fluctuation of nematic order parameter is 
not changed by random mass and random gauge field, 
which are both irrelevant at low energy. Therefore, it 
seems reasonable to expect an enhancement of dc ther- 
mal conductivity at low temperature if there are only 
these two kinds of disorders. However, when there is 
moderately strong random chemical potential, which is 
marginal, the extreme anisotropy fixed point is destroyed. 
Moreover, the nematic phase transition may become un- 
stable in the presence of such random chemical potential. 
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The extreme velocity anisotropy produced by the crit- 
ical nematic fluctuations may be probed by the heat 
transport measurements, since it leads to a remarkable 
enhancement of the low-temperature thermal conductiv- 
ity. Apart from transport measurements, such extreme 
anisotropy can also show its existence in angle-resolved 
photo-emission spectroscopy (ARPES) experiments. In- 
deed, the currently known value of the velocity ratio 
va/vf in d-wave cuprate superconductors was extracted 
from both heat transport [34[ and ARPES measurements 
35] . Unlike heat transport experiments that can only es- 
timate the ratio v\/vp, the ARPES measurements are 
able to determine the Fermi velocity vp and the gap 
velocity va separately 35]. From the solutions of RG 
equations, we know that, the extreme velocity anisotropy 
emerges because va is driven by the critical nematic fluc- 
tuations to drop rapidly down to zero at large I but vf 
is driven to decrease very slowly. It should be possible to 
detect the extreme anisotropy by measuring Dp and va 
separately by means of the ARPES experiments. 

In addition, the effects of the nematic order parameter 
fluctuations can also be reflected in the single-particle 
spectral function of nodal quasiparticles. Kim et al. [l8j 
investigated this issue and found two important features: 
strong angle- dependence of quasipaticle scattering, and 
an enhancement of velocity anisotropy. These predicted 
features of the fermion spectral function are expected to 
be tested by ARPES experiments. 

We next would like to remark on the disorder effects on 
the polarization function and the final results. In our cal- 
culations, the polarization function is obtained from the 
bubble diagram shown in Fig. [I] since including inter- 
nal disorder scattering line will introduce an additional 
suppressing factor gvp. To justify this approximation, 
we now make a qualitative analysis on the disorder ef- 
fects. It is well known that disorder scattering generates 
a fermion damping rate 70, which shifts the energy of 
Dirac fermions from lj to uj + 270- In the presence of 
finite 70, it seems possible to get a full analytical expres- 
sion for the polarization function Ii{q x , q y , e) only in the 
static (e = 0) limit. Following the computational proce- 
dures given in (3(| [37j , we have 



In order to simplify this expression, we now con- 
sider the low-energy regime, |q| < 70, which 

leads to arctan 1 



(70 1 \J /x ( l - x )( v fQx + V Aly)) 



7 1 y x(l — x){v 2 F q 2 + v\q 2 ). Using this approximation, 
the polarization function is found to be 



11(^,^,70) 



&-k 2 va!o 



ill +?«)• 



(72) 



Substituting this new polarization to the nematic prop- 
agator and then performing the same RG calculations, 
we found that the qualitative results presented in Sec. |V] 
do not change. We thus conclude that it is justified to 
neglect disorder effects in the polarization function for 
weak disorders. Admittedly, when disorders are strong 
enough to cause Anderson localization, the RG approach 
utilized in our manuscript is no longer applicable and a 
new RG scheme is needed to deal with the vertex correc- 
tions generated by disorder scattering. 

In this paper, we have considered only the coupling be- 
tween disorders and fermionic nodal quasiparticles since 
we are mainly interested in the disorder effects on the 
RG flow of fermion velocities. In practice, the effects 
of disorders on the nematic transition are more compli- 
cated. For instance, there might be quenched disorders 
that couple directly to the nematic order parameter. This 
issue was briefly discussed by Kim et al., who argued 
18] that quenched disorder may smear the symmetry- 
breaking type quantum phase transition thereby produc- 
ing a glassy state. It is currently unclear how the nodal 
quasiparticles are influenced by such kind of disorders. 
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n((7z,%,7o) 



1 



- 1 dx 2^/x(l~x)v F q 2 x 
2w 2 VFVA Jo . A,2 ~2 _u „2 ~2 



x arctan 

+{q x < — > %)■ 



"Hi + V A% 

(lo 1 ^(1 - x){v 2 F q 2 + v 2 A q 2 )) 

(71) 



Appendix 

In order to maintain a self-consistency of this paper, 
here we provide a detailed calculation of the fermion self- 
energy (16). Following Ref. [19j . we can define 



SWtf) = / 70^(Q + K)G{Q)K (^±r^) K (j^j , (73) 

where K = (k,w) and Q = (q, e) are 3-momenta. Here /C(y) is an arbitrary function with /C(0) = 1, and it falls off 
rapidly with y, e.g. /C(y) = e~ y . However, the results are independent of the particular choices of K,(y). It is easy to 



identify that, 



G(Q) = 
F(Q + K) = 



1 



n(q,e)' 

1 i(e + lo) - v F (q x + k x )r z + v A (q y + k y )r x 



N f (e + w) 2 + v 2 F {q x + k x f + v 2 A (q y + k y f 
Expanding F{Q + K)IC ( j^jp^ ) at Q + K = Q, and retaining the first order, we have 

•(q + k) 2 * 



F{Q + K)IC 



A 2 



K„ 



dF(Q + K)^ /(q + k) 



A 2 



dF(Q) 

Then the self-energy becomes 



+ F(Q + K) 



2 9 M y, ( q 



2(q + k)^, /"(q + k 



A 2 



\2\ 1 



A 2 



A 2 



A 2 



d 3 Q 
(2tt) 3 



which leads to 



dA 



d 3 Q 



(2tt) 3 



A 2 



-4q 2 0F(Q) 



2X 2o /a 2 

+ F{Q)G{Q)^1C I J 



A 2 



A 2 



4q 2 g M 
A 4 



4F(Q)^]G(Q)/C(^)/C'(^) 



F(Q )G (g)[/c(^)r'(^) +r 2 (^) 



It is convenient to introduce the following cylindrical coordinates, 

— yA(vFX,cos8,sm9), 

Qfi = ("fi, cos#, sin#), 

= 2/A(0,cos#,sin6>), 

q^ = (0,cos#,sin#), 

d 3 Q = y 2 A 3 v F dxdyd9. 



It is straightforward to obtain 



where 



F(Q) = 
G(Q) = 

g- 1 



1 / ix — cos9t z + (va/vf) sm0T x 
NfVF \ i 2 + cos 2 9 + (va/vf) 2 sin 2 9 

-L. = 1 6„ 4 e(x, e ), 



\jx 2 + cos 2 6* + {va/vf) 2 sin 2 
x 2 + sin 2 
a; 2 + sin 2 9 + (va/vf) 2 cos 2 9 



Since F and G are homogenous functions, 



HQ) = ~x F (Q)> °(Q) = ^x°(Q)- 
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We thus have 



dX nni (K) v F K^ 



dA 



(2tt) 3 



p2TC pOO 

dx I dO 



oo JO 



dy{ 



±y%F{Q)\ G{Q)JC(y 2 )K,'(y 2 ) 



-4y%F(Q)G(Q) [fC{y 2 )K,"{y 2 ) + K,'\y 2 )] } 
vfK^ 



8tt 2 



r 2ir 

dx I d6 



oo JO 



4 ^5^ - 4^F(Q)l G(Q) / ydy/C(y 2 )/C V) 



/■oo 

-4g„F(Q)G(Q) jf V 3 dy[K,(y 2 )JC"{y 2 )+K, l2 {y 2 )\). 



After integrating y out, we find 



f 

Jo 



y A dy lC(y 2 )JC"(y 2 )+)C /2 (y 2 ) 



ydylC{y 2 )lC'{y 2 ) = -. 



Therefore, 



A 



dA 



8tt 2 
vfK^ 



8tt 2 Jo dQ^ 

Formally, the fermion self-energy function can be expanded as 

d^nm(K) 



^ dx r dd{ + q^F(Q)] G(Q) - q»F(Q)G{Q)X 

J-oo Jo L L oQ„, 

f dx I 

J-oo Jo 



2 \e d -m G{ Q). 



dlnA 

When Kq = lu and Qq = vpx, we finally have 



Ci(-iw) + C 2 v F k x T z + C 3 v A k y T x . 



Ci(-iw) 



VfOJ 

8^ 



oo p2tt 

dx / d6 

-oo Jo 



dF(Q) 



G(Q) 



^r*! (do 

N f n 6 v F J_ oc J 
2(v A /v F ) 



8vfx 

- 71 i(x 2 + cos 2 9 + (v A /v F ) 2 sin 2 9) - 2x (ix - cos 9t z + (v A /v F ) sin 9t x ) 



(82) 



(83) 



(84) 



(85) 



(x 2 + cos 2 9 + (v A /v F ) 2 sin 2 9) 2 



iV/TT 3 

which directly leads to 



- ^ x 2 -cos 2 9-(v A /v F ) 2 sm 2 e . 

dx j d6— — ; ; , 2 , 2 g{x,0)(-lU}), 



o (x 2 +cos 2 + (v A /v F ) 2 sm 2 d) 



2(v A /v F ) 
Nfir 3 



oo p2ir 



dx / dO 



C2 and C3 can be obtained similarly. 



x 2 - cos 2 9 - (va/vf) sm # 
(x 2 + cos 2 9 + (v a /vf) 2 sin 2 #) 2 
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